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Effect of Electron Screening on the Collapsing Process of 
Core-Collapse Supernovae 
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ABSTRACT 

By using an average heavy nuclei model, the effects of the electron screening on electron 
capture (EC) in core-collapse supernovae are investigated. A one-dimension code based 
on the Wsl5M0 progenitor model is utiHzed to test the effects of electron screening 
during the collapsing process. The results show that, at high densities, the effects of 
EC on electron capture becomes significant. During the collapsing stage, the EC rate 
is decreased, the collapse timescale is prolonged and the leakage of the neutrino energy 
is increased. These effects lead to an appreciable decrease in the initial energy of the 
bounce shock wave. The effects of electron screeening in the other progenitor models 
are also discussed. 

Key words: nuclear reactions, nucleosynthesis, abundances - stars:evolution - super- 
novaeigeneral. 



1 INTRODUCTION 

Supernovae explosion is one of the most violent events in our 
universe. Their explosion mechanism is an old problem, but 
has not been understood completely ( see e.g. Woosley & 
Heger 2007; Martmez-Pinedo 2008 for reference) . Generally, 
a massive star(M ^ IOMq) proceeds through all burning 
' stage from Hydrogen to Silicon, finally leading to an iron 
, core in its center(Blanc & Greggio 2008). Electron capture 
(EC) causes the number density and the degenerate pres- 
sure of electron to decrease with a lot of neutrino energy 
loss, which lead to the accelerating collapse of the star till 
the central density reaches the maximum (two-three times 
the nuclear density). Later on, the infalling outer core col- 
lides with the stiff inner core, and then the bounce shock 
is produced in the vicinity of the boundary between the 
inner core and outer core. The initial temperature of shock 
wave is even higher than 10^^ K. At such a high temperature, 
photon energy is much larger than the binding energy of nu- 
cleus , and the iron nuclei are photodisintegrated into pro- 
tons, neutrons and electrons behind the shock wave(Janka 
2007). Unfortunately with the current numerical simula- 
tions, a self-consistent treatment of one-dimension super- 
novae model does not yet lead to successful explosion due to 
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the energy insufficiency, while two-dimension models show 
some promise (Ifix et al. 2003; Burrows et al. 2006). 

ft is well known that the weak interactions, espe- 
cially the EC and beta decay, are essential to the evolu- 
tion of supernovae. Many basic investigations at this aspect 
have been done by Bahcall(1962,1964), Fuller, Flower, New- 
mann (1982a; 1982b), Aufderheide et al. (1994), Langanke, 
Martinez-Pinedo (1999; 2000) and so on in the last decades. 
Fuller et al. (1982b) have ever mentioned the screening cor- 
rection but they did not make a detailed calculation. Later, 
Hix & Thielemann (1996) and Bravo & Garci'a-Senz (1999) 
considered the screening correction on the silicon burning 
and nuclear statistical equilibrium (NSE), and they found 
that screening effect is significant. Luo & Peng (1996) in- 
vestigated the effect of screening on EC in the supernovae 
environment by using Fuller et al. (1982a, b) method (i.e. 
so called shell model brink hypothesis). Their results show 
screening can reduce the EC rate by 10 — 20 per cent at 
high density. More recently improved EC rate and screening 
potential were adopted in Liu, Zhang & Luo (2007). The 
screening affects the EC mainly in three aspects. (1) The 
screening changes the Coulomb in- wave function of the elec- 
tron; however, it can be neglected because the screening po- 
tential is much less than the average energy of the electrons. 
(2) The screening reduces the energy of the electrons in the 
capture reaction. (3) The screening decreases the number 
of the high-energy electrons with energy greater than the 
threshold energy of EC. As a result, the reaction threshold 
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energy increases. Compared with screening on EC, there are 
at least two differences for beta decay. (1) The electron en- 
ergy of beta decay is decided by the energy difference be- 
tween parent nucleus and daughter nucleus (including rest 
mass), but the electron energy in EC can be much larger 
than that of the capture threshold. (2) In the inner core, 
beta decay is prohibited due to the Pauli exclusion princi- 
ple (inhibition degree has been given by Peterson & Bahcall 
1963); at the outer core, beta decay gets permission, but the 
evolution of inner core is quicker than that of outer core at 
core collapsing stage (which is less than 0.5s). So we here 
neglect the screening effect on the beta decay. Langanke et 
al. have suggested the screening should be considered in the 
simulations (Langanke & Martmez-Pinedo 2003), but up to 
now the screening effect on the total explosion process has 
not been done. Early numerical simulations show that the 
onset of the SN explosion depends on the pre-supernovae 
model and the evolution mode that is especially sensitive 
to the physical parameters input, which is closely related to 
the electron fraction and the weak interaction rates (Heger, 
Woosley & Martmez-Pinedo 2001). Therefore, it is impera- 
tive to obtain the EC rates with high precision. As an actual 
physical input, screening should not be ignored since a se- 
ries of important parameters, such as neutrino energy loss, 
collapsing timescale etc., is changed. In this paper, we in- 
vestigate the detailed effects of electron screening on EC in 
core- collapse supernovae. 



2 PRE-SUPERNOVAE MODEL AND 
COMPUTATIONAL APPROACH 

We perform numerical simulations by using a modified ver- 
sion of the one-dimension code developed by Y.R. Wang, 
S.C. Zhang, W.Z. Wang, Z.H. Xie (WZWX 2003) in 1996 
. In this code the general relativistic hydrodynamic equa- 
tions are adopted (May & White 1966); the hydrodynamics 
method is smoothed particle hydrodynamics (SPH) referring 
to Benz (1991); equations of state are similar to Lattimer 
& Swesty (1991) and Cooperstein & Wambach (1984); the 
method of neutrino transport is provided by Suzuki (1994). 
The pre-supernovae model we choose is the WslSAf© model 
with an iron core of 1.38M0 (Woosley & Weaver 1995). The 
grid domain includes 1.6 Mq, which is divided into 96 mass 
layers. Since most of nuclei involved in supernovae environ- 
ment are on the unstable and excited states, it is difficult 
to get a precise description of each energy level, especially 
for heavy nucleus whose excited states are almost continu- 
ous. For a simple consideration, we assumed that the matter 
consists of four typical particles including free protons, free 
neutrons, a-particles, and heavy nuclei, i.e. the so called 
'four-particle model'. Such four types of particles can well 
represent the whole property of pre-supernovae (Lattimer 
& Swesty 1991; Arcones Janka & Scheck 2008). Similar to 
the other authors, detailed numerical model indicates the 
shock is unable to rush out of the iron core because of too 
much energy loss in the iron photodisintegration. So here 
we mainly investigate the screening effect in the collapsing 
process of supernovae explosion. 

Usually, the capture rates for the nucleus {Z, A) in ther- 
mal equilibrium at temperature T is given by a sum over the 



initial parent states and the final daughter states(Pruet & 
Fuller 2003), 
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where Ji and Ei are the spin and excitation energy of the 
parent states, respectively, fes is the Boltzmann constant 
and G[Z, A, T) the nuclear partition function. The /t-values 
are related to GT and Fermi transition matrix elements, fif 
is the phase space integral of electron. But for the 'four- 
particle model' we adopt, matter is composed of protons, 
free neutrons, a-particles, and heavy nuclei under NSE. Only 
proton and heavy nuclei can capture electrons, so the total 
EC rate A = Ap -I- Ah, where \p and \h are the EC rates for 
proton and heavy nuclei, respectively. Considering the core 
of the pre-supernovae is composed of iron element, the dom- 
inant influence of electron screening is decided by the heavy 
nuclei. Precise calculation of the EC rate is usually based 
on the nuclear shell model, i.e formula ([1]), but the short- 
age of energy level quality at high temperature and diverse 
electron rates for so many different nucleus will bring com- 
plication and some uncertainty. Therefore, here the Fermi 
gas model is adopted to describe the average property of 
the heavy nuclei. The EC rate is defined by A = —{dne/dt), 
when the number density of electron decreases A > 0; oth- 
erwise A ^ 0. According to Bethe et al. (1979), the EC rate 
for average heavy nuclei is expressed 
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where = pNAYe,np = xxhpNa,Pp = {2mpfip 
the density, Na is the Avogardro constant, fip is the chemi- 
cal potential for proton. Ye is the electron fraction, nip the 
mass of proton, x = Z/A, Z, A are the mass number and 
charge of average heavy nuclei, respectively, xh is the frac- 
tion of the average heavy nuclei, Se and Sp is the energy 
of electron and proton, respectively, /e = {1 + exp[(ec — 
fie ~ meC^) / {kBT)]}~^ is the Fermi-Dirac distribution func- 
tion for electrons. Because of the conservation of energy, 
Ee + £p = £n + £v, equation ([SJ becomes 
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where Q is the reaction threshold energy. Q — \Qif\ if 
Qif < -0.511MeV, otherwise Q = O.SUMeV. Q,/ = 
{MpC^ - MdC^ + - Ef), Mp and Ma are the mass of 
parent and daughter nucleus, respectively. Ei and Ef the 
excited energy for parent and daughter nuclei, respectively. 
Note that here parent and daughter nuclei are proton and 
neutron in bound state of nucleus, so Qi/ = /i„ 4- A„ — Sp, 
which can also assurance the energy of neutrino is not less 
than zero. Proton energy Sp = fip — Ap, therefore 



Hp - ifJ-n + A„ + £v - fie) 
fie - P- - An ~ Ev ^ A - , 



(4) 



where p, = p„ — pp, A — pe — p — A„ — 
1.15p'[q''^ y^JiT/Yeidp/dpey^* is the maximum energy of 
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emissing neutrino, average energy of neutrino ~ (3/5)A, 
so Ap ^ (2/5) A. A„ = mm(3, max(0, A/2)), where 3MeV 
is a refereed value given by Bethe et al. (1979). Some values 
at initial moment are listed in Table 1. 

In the high density gas, screening electron cloudy is 
formed. Screening reduces the energy of the electron and en- 
hances the threshold energy, so the EC rate will decrease. We 
adopt the method similar to Liu et al. (2007), i.e., we ignore 
the effect screening changes the Coulomb in-wave function 
of the electron, only consider the effect on the electron en- 
ergy and threshold energy. So the capture rate in the strong 
screening is rewritten as 



Ah 
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where D is the modification to the threshold energy, Vs is 
the screening potential(Fuller et al. 1982b, Itoh et al. 2002): 



D = 2.94 X W-^Z^^^ipYeY'^iMeV) 



= 7.525 X lO"''^ 



i^^pej J{MeV) 
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E'^^o'^^J^'""' = 0.5(logr, + 3),u = l/25(i? - 25), r» = 
1.388 X 10-^(A/^p6)^''^ 



(2002), equation dH) is valid for 10"^ ^ ^ 10"^ (i.e. 
p ^ lO^^g cm"'^ ), which is usually fulfilled in the super- 
novae environment. 

For the EC on proton, the most probable energy of in- 
teraction _Bo — Ef- Because of i5o >> Ep and kT > Ep, 
where Ep is the Coulomb energy for proton, strong screening 
is invalid and weak screening is valid (Salpeter & von Horn 
1969). According the method of Kippenhahn & Weigert 
(1990), Bahcall, Chen, & Kamionkowski (1998) and Kippen- 
hahn & Weigert (1990), weak screening factor is exp(a;ovr7;) 
for XQ « 1, where xo = r^jro ^ iZiZ2/200Eo){Cg/T7y^^ 
(note that here Eq is in KeV); r^, td are the classical 
turning-point radius and Debye radius, respectively; C, = 
yJY^.{XiZ'^ /Ai + XiZi/Ai), Xi,Zi, and Ai are the mass 
fraction, charge, and mass number, respectively; Ty is tem- 
perature in unit of lO'^K, = {m/2Y^'^[{ZtZ2e^) / [nEy^)], 
where m is the reduced mass. With rough estimate for 
the core region of supernova, T ~ l(fK, C, ~ 10, and 
p = 10® ~ lO^^gcm"^, weak screening correction is not 
more than 0.001. On the other hand, proton fraction is much 
smaller than heavy nuclei, so screening correction for free 
proton capture is not important comparing with that on 
heavy nuclei. 



3 SIMULATION RESULTS 

Because there are many output parameters in the simula- 
tion, in Table 1 we only tabulate some important parameters 
(only 10 layers are listed) at the initial moment of simula- 
tion (about 0.27s before bounce). One can see that the mass 
number and charge of average heavy nuclei, density, elec- 
tron fraction, temperature in the pre-supernovae model and 




Figure 1. Neutrino energy loss rates at 0.27s before bounce. Solid 
curve is neutrino energy loss rate without screening, while dashed 
curve is that with screening in the different layers. 




Figure 2. The electron capture rate at different collapsing time. 
Solid curves are EC rates without screening, while dashed curves 
are that with screening in the different layers. The symbols 1, 2, 
3, 4, 5, 6, 7,8, 9 represent the time of -250.0, -200.0, -100.0, -50.0, 
-25.0, -5.0, -1.5, and 6ms for the case of without screening ( 
'— ' means the time is before bounce, ms is the abbreviation for 
milliseconds). Squares denote the neutrino trapping positions at 
the corresponding time. 



screening potentials at the different layers. Here, A and A' 
represent the EC rate without/with screening, respectively. 
It is shown in Table 1 that the mass number and charge of 
heavy nuclei decrease from the center to the outer layer. The 
screening potential is mainly dependent on the density, so 
both D and Vs decrease monotonically with densities. A is 
not only a function of Z, A, Ye, and T, but also a function of 
Ap,Aif,Xp and xh, so it is not monotonic to the mass layers. 
With screening effect into consideration, it is easy to find 
that A' is always smaller than A. 

The neutrino leakage and diffusion process are the key 
factors on the explosion energy, and the energy loss of neu- 
trino exceeds 90 per cent of the total in the whole evolution 
process of the star. Neutrino energy loss rate is an important 
and complicated parameter to affect the explosion energy 
and evolution, which is related with the weak interaction 
rate, neutrino energy, and transport equation and so on. 
The change of EC rate must influence the neutrino energy 
loss rate. Fig. [T] shows the comparison of neutrino loss rate 
with/without screening at 0.27s before bounce. Solid curve 
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Table 1. Some important parameters at the initial moment0.27s before bounce. 



j A Z p Ye T D Vs fip Pn Ap A„ A A' 

gcm,-^ K MeV MeV MeV MeV MeV MeV s-^cm.-^ s-'^cm,-^ 



1 


61.6 


26.1 


8.3EH 


-09 


0.42 


8.1E+09 


0.39 


0.14 


-12.54 


-6.72 


0.37 


0.92 


6.28EH 


-30 


5.41E-I-30 


10 


60.1 


25.8 


5.3EH 


-09 


0.43 


7.9E+09 


0.33 


0.12 


-12.04 


-7.06 


0.32 


0.79 


1.86EH 


-30 


1.65E-I-30 


20 


58.6 


25.4 


3.4EH 


-09 


0.43 


7.7E+09 


0.29 


0.10 


-11.47 


-7.46 


0.33 


0.82 


1.13EH 


-30 


1.02E-I-30 


30 


57.1 


25.1 


2.2EH 


-09 


0.44 


7.5E+09 


0.25 


0.08 


-10.83 


-7.93 


0.39 


0.98 


1.14EH 


-30 


1.05E-(-30 


40 


55.0 


24.6 


1.3EH 


-09 


0.45 


7.3E+09 


0.21 


0.07 


-9.95 


-8.60 


0.54 


1.35 


1.71EH 


-30 


1.59E-I-30 


50 


52.9 


24.1 


7.5EH 


-08 


0.46 


6.8E+09 


0.17 


0.06 


-9.18 


-9.18 


0.66 


1.65 


1.50EH 


-30 


1.40E-I-30 


60 


52.1 


23.9 


3.4EH 


-08 


0.46 


7.1E+09 


0.13 


0.04 


-9.19 


-9.19 


0.47 


1.17 


1.80EH 


-29 


1.73E-I-29 


70 


49.6 


23.2 


1.2EH 


-08 


0.47 


6.1E+09 


0.09 


0.03 


-9.13 


-9.13 


0.31 


0.77 


7.20EH 


-27 


6.90E-I-27 


80 


43.8 


21.8 


2.7EH 


-07 


0.50 


4.7E+09 


0.06 


0.02 


-8.89 


-8.89 


0.17 


0.41 


9.08EH 


-25 


8.77E-I-25 


90 


43.3 


21.3 


2.0EH 


-06 


0.50 


2.9E+09 


0.02 


0.01 


-8.80 


-8.80 


0.04 


0.11 


3.57EH 


-22 


3.48E-I-22 



note, j denotes the layer number. 
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Figure 3. The explosion energy with/without screening versus 
collapsing time. The notes are the same as those in Fig.l. 




R(cm) 



Figure 4. The collapsing velocity versus radius at bounce. Shock 
is produced at sound speed position, and radius of PNS appears 
at the vicinity of maximal collapsing velocity at the bounce. The 
notes are the same as those in Fig.l. 



denotes loss rate without screening, dotted curve denotes 
that with screening. We find neutrino energy loss rate de- 
creases generally, at some region the modification is more 
than 5 per cent (It closely depends on EC rate), but for 
outer part of iron core, it has hardly changed. 

Fig. [2] shows the evolution of EC rate with/without 
screening at different moment during the collapsing stage. 
One can find (i) the EC rate with screening is always smaller 



than that without screening, which is caused by the decrease 
of the EC rate and the delay of collapse time-scale. At the 
initial stage (such as symbol 1), the difference is mainly 
caused by the different EC rates, and at the later stage (such 
as symbols 5,6,7), the difference comes from a different col- 
lapse history, while at bounce moment, for the case without 
screening (symbol 8), the EC rates are similar to that at 
moment 6ms later for the case with screening, (ii) In the 
inner part of the core, both the EC rates and the differ- 
ence between the case with screening and the case without 
screening increase rapidly with time. The reason is that the 
densities increase quickly with the collapsing process. The 
higher the density, the larger the screening potential. How- 
ever the screening effect is decided not only by the potential, 
but also by the Fermi energy of electrons. When the den- 
sity is large enough, the Fermi energy becomes much larger 
than the potential energy. Correspondingly the screening ef- 
fect will decrease. (iii) In the outer part of the core, the EC 
rates change comparatively slow. At the edge of iron core 
(1.38Mq), there is an infiexion because both the fraction of 
heavy nuclei and the densities break obviously. At the outer 
envelop, EC rate is almost stable since the time is just 0.3s. 
(iv) Note the squares symbols in Fig. 2, from left to right, 
they denote neutrino trapping critical position at -5.0, -1.5 
and 0.0ms before bounce, respectively. We can find the high 
density region extents rapidly with time. Because of the in- 
verse process of EC rate [ve -f (A, Z) ^ + (A, Z -f 1)] 
enhances with density, above critical density {ptrapping is 
about 3 X lO^^pcm"^), an equilibrium of neutrino will be 
reached and the lepton fraction (including electron fraction 
and neutrino fraction) will keep as a constant. Surely, at 
such situation our method for screening modification is in- 
valid to change of lepton fraction, but screening is still exist 
and keep effect in the reversible reaction. 

We also compare the collapsing time and the radius 
of protoneutron star (PNS) with screening/without screen- 
ing(as shown in Fig. 3 and 4). We find (i) the collaps- 
ing time is prolonged; t=0.267s without screening while 
t=0.272s with screening. The reason is that screening de- 
creases the capture rate. This makes the electron number 
density and degenerate pressure to drop a little more slowly 
than the fiducial case. During the collapsing stage, the to- 
tal pressure is dominated by the electron degenerate pres- 
sure, so that the decrease of total pressure also become more 



Effect of Electron Screening on CCSNe 5 



slowly and the collapsing velocity decreases comparatively, 
(ii) The initial energy of shock wave decreases. The initial 
energy of shock wave at about 0.8ms after bounce is 1.06 
foe (lfoe=l X \(f^erg) without screening and l.Olfoe with 
screening. Their difference is 0.05 foe, which is 5 per cent 
of the total. On the one hand, it could be that the calcula- 
tion with screening results in a different radial profile of the 
star at bounce when compared with a calculation without 
screening. The different profile could result in differences in 
the propagation of the shock wave. Figure 4 shows that the 
velocity of the collapsing matter is different for radii above 
the shock position; on the other hand, it could be the total 
neutrino energy loss increase due to the time delay. The to- 
tal neutrino energy loss Q-ut = QuAMjAt , where 
t is the time, At is the time step, AAf,- is the mass of the 
jth layer. Despite the screening effect decreases the capture 
rate and the neutrino energy loss rate, the total energy loss 
is also related to the collapsing timo.(iii) Radius of PNS is 
hardly changed, but the radius of outer part of iron core 
increase slightly. It is favorable to the success of supernovae 
explosion. 



4 CONCLUSION 

In this paper, we have considered screening effect on EC, 
the corresponding effects on the collapsing process, and on 
the initial shock wave energy of core-collapse supernovae. 
Our results show that the screening effect is appreciable. 
Since there are both advantages and disadvantages to the 
final explosion, more detailed simulation is needed by us- 
ing more concrete equation of state and progenitor model 
(Hix et al. 2003). In order to investigate the exact effect of 
screening, the concrete nuclide and the other methods are 
needed. In recent yeaxs, many methods such as the large- 
scale shell model and the pn-random phase approach are 
widely investigated (Nabi & Sajjad 2008). Improved results 
are universally one order of magnitude smaller than Fuller 
et al. (1982a,b) results. If the improved EC rates arc con- 
sidered, the screening effect should increase comparatively. 
Furthermore, for two or three dimension numerical simu- 
lations of core-collapse supernovae, our method for electron 
screening is also valid. The explosion mechanism of the core- 
collapso supernovae has been investigated extensively in the 
last four decades and significant progress hcis been obtained, 
but some of the most fundamental questions are still unan- 
swered. It is necessary to include the effect of screening in 
later detailed simulations. 
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